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In a coalescence plus independent fragmentation approach we calculate the pr spectra of the main 
hadrons: tv, K,p,p, A in a wide range of transverse momentum from low pr up to about 10 GeV. 
The approach in its main features was developed several years ago at RHIC energy. Augmenting 
the model with the inclusion of some more main resonance decays, we show that the approach 
correctly predicts the evolution of the pr spectra from RHIC to LHC energy and in particular 
the baryon-to-meson ratios p/tTjP/tt, A/A that reach a value of the order of unit at pt ~ 3GeV. 
This is achieved without any change of the coalescence parameters. The more recent availability of 
experimental data up to pr ~ 10 GeV for A spectrum as well as for p/n and A/K shows some lack 
of yield in a limited pr range around 6 GeV. This indicates that the baryons pt spectra from AKK 
fragmentation functions are too flat at pr ^ 8 GeV. 

PACS numbers: Valid PACS appear here 


I. INTRODUCTION 

Heavy-Ion collisions at ultra-relativistic energies gener¬ 
ates the creation of a new state of matter with a temper¬ 
ature T significantly above the cross-over temperature 
Tc ~ IGOMeV expected for the transition of hadronic 
matter to the Quark-Gluon Plasma (QGP) [I|. In the last 
decade a consensus seems to have emerged that quarks 
and gluons are indeed deconfined for a short amount of 
time in the fireball created, and that this quark gluon 
plasma behaves like a very good liquid with small vis¬ 
cosity over entropy density ratio 77 /s [2~ki . One of the 
most surprising observations at RHIC |52l| have been the 
strong enhancement of the baryon over meson ratio in a 
wide range of 3 GeV ^ Pt ^ 6 GeV and a significantly 
larger (about 40 — 50% ) ellij^ic flow V 2 of baryons with 
respect to that of mesons [sHtI, at variance with expec¬ 
tation from both hydrodynamics and jet-quenching frag¬ 
mentation. However soon after these first data it was 
realized that an hadronization process through quark co¬ 
alescence can naturally explain both features of hadron 
spectra appearing in the same pr region [Ml- 

The idea of quark coalescence as hadronization from a 
QGP where the quark (anti-quark) needed to hadronize 
are already present and there is no need to create them 
from vacuum, was initially suggested in Ref.s [S in 
with the approach known as ALGOR model. However, 
there the main focus was on particle yields and ratios 
and not the different kinematics between baryons and 
mesons that explains their pattern for pp spectra and 
elliptic flow. 

After the first papers applying the coalescence mech¬ 
anism for hadron production at RHIG there 

have been mainly two further studies.In one case the co¬ 
alescence approach has been extended to include finite 
width that accounts for off-shell effect which allows to in¬ 
clude the constraint of energy conservation [I^i • This 


is seen not to change the general feature of the coales¬ 
cence mechanism, the baryon/meson enhancement and 
approximate quark number scaling of V 2 , but shows that 
at low Pt, when energy conservation plays a stronger 
role, a scaling with rap = {pp + is favored, at 

least for direct hadron production. More recently, an 
approach that from one hand enforce algebraically the 
unitarity similarly to ALGOR approach, but for the mo¬ 
mentum distribution use essentially the prescription of 
Ref.s d, [ 1 ^ has been developed in [2l| . The price 
paid to have both aspects is the need to add a global con¬ 
stant in the coalescence probability in such a way that 
all the quarks recombine and moreover the quark distri¬ 
bution functions are fitted by mean of six parameters to 
describe hadron spectra. Nonetheless the result are simi¬ 
lar to @,[13 except a better description of the low pp re¬ 
gion, where however the present coalescence scheme can 
be anyway questionable and the success of the descrip¬ 
tion of the spectra down to low pp may simply come from 
the large flexibility of the fitted parton distribution. 

In this paper, we employ the coalescence approach de¬ 
veloped in [io, [ 2 ^ that is based on the same basic idea of 
phase-space quark coalescence as in but being solved 
by a Monte Garlo approach allows to include a 3D geome¬ 
try as well as radial flow correlation in the partonic spec¬ 
tra and the effect of the main resonance decays. The last 
is mainly discussed in [13 and allows to have a reasonable 
description of the spectra also at low momenta that are 
usually dominated by feed-down from resonance decays. 
However the model does not implement a unitarity con¬ 
servation that guarantees hadronization of the full bulk 
of particles, also energy conservation and the finite width 
(off-shell) effect are not included as in BO which at 
low momentum play a role. Therefore the model has to 
be considered really applicable only at pp > 1.5 GeV, 
even if the description of the spectra or baryon/meson 
ratios remains approximately good also at lower pp, es¬ 
pecially for pions that are dominated by decay feed-down 
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at low pt- In fact the present dynamical formulation of 
the coalescence used remains fully meaningful when the 
single process of hadronization is small and hence it is 
not justified where the bulk of the particles resides. 

We report in this paper the results applying the ord¬ 
inal model of coalescence plus fragmentation in |10l | 
at both LHC Pb+Pb ^/sWn = 2.76 ATeV and RHIC 
Au+Au ^/sNN = 200AGeV energies. This last case 
was partially present also in , but some more res¬ 
onance decays have been added, the update Albino- 
Kniehl-Kramer (AKK) fragmentation function and also 
the A momentum spectra have been evaluated. Further¬ 
more the availability of experimental data allows for a 
comparison in a wider pT range. We see that the same 
model is able to correctly predict the spectra of the main 
hadrons (tt, K,p,p, A), and correctly account for the evo¬ 
lution of such spectra from RHIC to LHC energy and 
in particular the pr dependence of the pjTr,pfn, Aj 
ratios. 

The paper is organized as follows: in Sec. H we de¬ 
scribe the general formalism of a covariant coalescence 
model for mesons and baryons and how it is numerically 
solved by means of a Monte Carlo method. How mini¬ 
jets and the quark-gluon plasma partons are determined 
is described in Sec. HI. Results for the transverse mo¬ 
mentum hadronic spectra obtained from the coalescence 
model are given in Sec. IV for RHIC energy and in Sec. 
V for LHC energy. Finally, we conclude in Sec. VI with 
a summary of present work and an outlook about future 
developments and applications of the parton coalescence 
model. 


II. COALESCENCE PLUS FRAGMENTATION 
MODEL 

The hadronization mechanism is well known to belong 
to the non-perturbative domain of the QCD dynamics 
and certainly a first-principle description of hadron for¬ 
mation has yet to be obtained. One of the most common 
approach dealing with hadronization, routinely used in 
nuclear and particle physics for inclusive hadron produc¬ 
tion, is the independent fragmentation in which a single 
colored parton a has to hadronize into the hadron h. For 
this purpose, fragmentation or parton decay functions 
Da{z) —>■ h have been defined and give the probability of 
finding hadron h in parton a with a momentum fraction 
z, 0 < z < 1. The cross section for inclusive hadron pro¬ 
duction can then be written as ah — era's) Da^h which is 
a convolution of the production cross section Ua for par- 
ton a with the fragmentation function Da^h{z)- Frag¬ 
mentation functions are not reliably calculable from first 
principles in QCD. However, they are observables and 
can be measured experimentally. 

Physically, the fragmentation of a single parton hap¬ 
pens through the creation of qq pairs, which subsequently 
arrange into color singlets forming hadrons. Such a 
scheme is based on the concept of QCD factorization. 


which separates the long- from the short-distance dy¬ 
namics. Notice that in a bulk QGP matter created in 
heavy-ion collisions the hadronization can occur on av¬ 
erage at a time scale th ^ ^ — lOfm/c (depending on its 
initial temperature) this means that the partons do not 
need a creation of qq at the moment of hadronization, 
but they can rescatter and generate a thermal medium 
expanding with a radial collective flow. At th partons 
have a certain abundance in phase space such that there 
is no need for the creation of additional partons through 
splitting or string breaking. The most naive expecta¬ 
tion for such a scenario is a simple recombination of the 
deconfined partons into hadrons. Indeed, there is exper¬ 
imental evidence that this is the correct description of 
hadronization, even long before a thermal occupation of 
parton phase space is reached (23l - l^ . 


A. Coalescence 


The approaches developed in S ES El, HI for the 
coalescence is based on the Wigner formalism originally 
developed for nucleons [2^. The spectrum of hadrons 
formed from the coalescence of quarks can be written as: 


(PNh f TT- (Ppi q f ( \ 

-Iff = s-j 11 

Ept.) (1) 


where dai denotes an element of a space-like hyper¬ 
surface, gn is the statistical factor to form a colorless 
hadron from quark and antiquark with spin 1/2. fq. 
are the quark (anti-quark) distribution in phase space. 
Jh is the Wigner function and describes the spatial and 
momentum distribution of quarks in a hadron. Eq. [1] 
for i = 2 describes meson formation, and for i = 3 the 
baryon one. The above formula can also be used for an- 
tibaryons by replacing quark momentum spectra by the 
momentum spectra of antiquarks 

The statistical factor gn takes into account the internal 
quantum numbers in forming a colorless meson from spin- 
1/2 colored quark and antiquark. For mesons considered 
here, i.e., tt, p, a;,iF, and AT*, the statistical factors are 
diT = dK = 1/36 and gp = gx* = 1/12. For baryons and 
antibaryons considered in present study, i.e. p,A,p and 
A, A, the statistical factors are gp = gp = gA = 1/108 
and 5A = ffA = 1/54. The coalescence probability func¬ 
tion fMixi,X 2 ;pi,P 2 ) ifHixi,X2,X3;pi,p2,P3))depends 
in principle on the overlap of the quark and anti-quark 
distribution functions with the wave function of the me¬ 
son (baryon) as well as the interactions of emitted vir¬ 
tual partons, which are needed for balancing the energy 
and momentum, with the partonic matter. Neglecting 
the off-shell effects the coalescence probability function 
is then simply the covariant hadron Wigner distribution 
function. 
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In the Greco-Ko-Levai (GKL) in Ref. 0 approach for 
a meson the Wigner function is a sphere of radius A^; 
in coordinate space and Ap in momentum space, these 
two parameters are related by the uncertainty principle 
A,,Ap = 1, 

9 TT 

fM{xi,X 2 \Pl,P 2 ) = -Y 

X e { Al-pl + Aml ^) 

where we have defined the quadri-vectors for the relative 
coordinates Xri = Xi — X2, Pri = Pi — P2 and the scalar 
Ami 2 = mi—m 2 . For a baryon we have a similar Wigner 
function expressed in term of appropriate relative coor¬ 
dinates: 

fB = - IxJ.) X e (a; - ip;,) 

0 (A^ — X^ 2 ) ^ ® i^p ~ Pr 2 ~ ^^ 1123 ) ) (2) 

and we have defined Xr 2 = {^ci + X 2 — 2 x 3 )/'/E and 
Pr2 = {pi +P2- 2 p 3 )jV 6 and Ami23 = mi +TO 2 - 27713 . 
We note that the Wigner function for the meson has only 
one parameter Ap = A~^ that we choose to have a mean 
square radius of mesons of about 0.8 fm which corre¬ 
sponds to a = 0.19 GeV. For baryons we have 

also one parameter, even if one could consider two differ¬ 
ent ApS for the two relative coordinates, for qqq baryons 
we have a = 0.33 GeV and for qqs baryons 

Ap = 0.38 GeV which corresponds to a slighter smaller 
radius for A as one can expect from the wave function of 
an harmonic oscillator that scales with the inverse square 
root of the reduced mass of the system. However we have 
checked that small variation of Ap does not strongly af¬ 
fect the slope of the spectra coming from coalescence. 

A specific feature of the GKL approach has been the 
inclusion of resonance decays to the pion, proton, K spec¬ 
tra. This was a way to include a minimal effect that cer¬ 
tainly is present in the final spectra observed. In this 
paper we have also included some more resonance with 
respect to the original work [l^ this however does not 
change significantly the results for the pt spectra espe¬ 
cially at Pt ^ 1.5 GeV. This can be envisaged because 
the resonance decay products mainly feed-down the low 
Pt region. The residual effect at intermediate pt can be 
reabsorbed by a rescaling of the Wigner function hadron 
parameter Ap by about a 10 — 15%. For the resonances 
the coalescence probability is augmented with a suppres¬ 
sion factor that takes into account for the Boltzmann 
probability to populate an excited state of energy AE at 
a temperature T, namely accordingly to the statistical 
model factor exp{—AE/T), with AE = Eh* — Eh and 
Eh* = (pt + and m*pj the mass of the resonance; 

of course also the degeneracy factors that come from the 
different values of isospin and total angular momentum 
are taken in account. 


B. Fragmentation 

In Ref. [ 13 , M El it has been clarified that at in¬ 
creasing Pt the probability to coalescence decreases and 
eventually the standard independent fragmentation takes 
over. It is therefore necessary to include also the contri¬ 
bution from the mini-jet fragmentation. This is done 
by employing the same parton distribution function that 
at high Pt > Po ^ 3 GeV are those that can be calcu¬ 
lated in next-to-leading order (NLO) in a pQGD scheme. 
However in AA collisions one must include also the mod¬ 
ification due to the jet quenching mechanism [H El. 
The hadron momentum spectra from the minijet parton 
spectra is given by: 

dAfjiad _ f j dNjet ^had/jeti^iQ ) 

d'^PTdy ~ d^PTdy U 


where z = phad/pjet is the fraction of minijet mo¬ 
mentum carried by the hadron and = {phad/2,z)‘^ is 
the momentum scale for the fragmentation process. For 
Dhnd/iefiz, Q^) we employ the AKK fragmentation func¬ 
tion |3ll . l32j| . The fragmentation is applied for partons at 
Pt > 3 GeV. It is also worth to mention that the inclu¬ 
sion of both coalescence and fragmentation does not lead 
to a double counting of hadron formed because a hadron 
at p^ from coalescence comes from partons at pT — Pt 
with n = 2,3, while from fragmentation the parton cre¬ 
ating it comes from a pt — 1.5 — 2p^. Therefore, for 
example considering the fragmentation for a hadron at 
Py ~ 3 GeV the quark that creates it has on average a 
Pt ~ 5 GeV but the probability of coalescence for such 
a parton is very small and it would produce hadrons in 
a very different region of p!^ ^ 10 — 15 GeV where the 
coalescence is certainly negligible. 

As in El the multi-dimensional integrals in the coales¬ 
cence formula as well as the fragmentation are evaluated 
by the Monte-Garlo method via test particles. Specifi¬ 
cally, we introduce a large number of test partons with 
uniform momentum distribution in phase-space. To take 
into account the large difference between numbers of 
thermal and minijet partons, a test parton with momen¬ 
tum Pt is given a probability that is proportional to the 
parton momentum distribution, e.g., dNq/d'^pT, with the 
proportional constant determined by requiring that the 
sum of all parton probabilities is equal to the parton 
number. With test partons, the coalescence formulas, 
for mesons and baryons can be re-written as 


9M Pq{i)Pq{j)d^^\pT - P*T - Pjt) 

fM{xi,Xj;p„Pj). (4) 


dNjH _ 
d^PT 

X 
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and 

93 ^ Pq{l)Pqij)Pq{k) 

- PiT - PjT - Pfcx) 
fB{Xi,Xj,Xk-,Pi,Pj,Pk)- (5) 

In the above, Pq{i) and Pq{j) are probabilities carried by 
ith test quark and jth test antiquark. 

The Monte-Carlo method allows us to treat the coa¬ 
lescence of low momentum partons on the same footing 
as that of high momentum ones. 


dNs _ 
fPv>T 

X 

X 


III. FIREBALL AND BARTON 
DISTRIBUTIONS 

The coalescence approach as developed till now to de¬ 
scribe hadron production in uRHIC’s is based on a fire¬ 
ball where the bulk of particles is a thermalized system 
of gluons and u, d, s quarks and anti-quarks at the tem¬ 
perature Tc = 165 MeV [ij which is about the tempera¬ 
ture for the cross-over transition in realistic lattice QCD 
calculation [s^. At high momenta, pr > 2.5 GeV, the 
distribution is taken to be the partonic spectra that un¬ 
dergone the in-medium jet quenching [33| . The longitu¬ 
dinal momentum distribution is assumed to be boost- 
invariant, i.e., a uniform rapidity distribution in the 
range y G (—0.5, -1-0.5). While these features are common 
to both RHIC and LHC, with the beam energy there is a 
change in the radial flow and self-consistently in the vol¬ 
ume of the bulk of the QGP, as well as at high pr there 
is a different initial distribution and jet quenching. 

In this paper, we consider Au-fAu at -^/snn = 200 
GeV at RHIC and Pb-I-Pb collision at ^snn = 2.76 
TeV at the Large Hadron Collider (LHC) for the central 
collisions (0 — 10%). To take into account for the quark- 
gluon plasma collective flow, we assume for the partons 
a velocity linear radial profile as Pt = Pmax^i where R 
is the transverse radius of the fireball. The quark and 
antiquark distribution up to 2 GeV are hence given by a 
thermal distribution with flow: 

dNq, q _ gq^qTrriT f jT{mT -pr ■ PtT Pq) \ 

dPrr d^pT ~ (27r)3 V T ) 

( 6 ) 

here gq = Pq = G indicates the spin-color degeneracy of 
light quarks and antiquarks, and the minus and plus signs 
are for quarks and antiquarks, respectively. The temper¬ 
ature is set to T = 165 MeV. Eq. m also applies to 
gluons after replacing pq by the gluon spin-color degener¬ 
acy pg = 16 and dropping the chemical potential. For the 
gluon mass, we take it to be similar to that of light quarks 
in order to take into account non-perturbative effects in 
the quark-gluon plasma. Because rUg < 2mq, we convert 
pp —>■ qq. Positions of partons in the transverse direction 
are taken to have a uniform distribution. Their longitu¬ 
dinal positions are then determined hy z = rsinhy, as 
we have assumed Bjorken correlation p = y. 


The radial flow Pmax and the volume V = ttR^t (in 
one unit of rapidity) could be in general considered as 
parameters evaluated accordingly to the typical value of 
lifetime of the QGP and, assuming a a constant accel¬ 
eration, we connect the radial expansion with the ra¬ 
dial flow i?_L = i?Q -I- O.bPrnaxT. So the radial flow and 
the volume are constrained imposing the total multiplic¬ 
ity dN/dy and the total transverse energy dEr/dy to 
be equal to the experimental data. The charged parti¬ 
cles multiplicity per unit of rapidity is dNch/dy ~ 1600 
and the transverse energy dEr/dy = 2100 GeV at LHC, 
dNch/dy ~ 670 and dEr/dy = 760 GeV at RHIC. This 
leads for a radial uniform expansion to Pmax = 0.37 and 
R± = 8.7 /to, t = 4.5 fm/c at RHIC , and Pmax = 0.60 
and R± = 10.2/to, r = 7.8 fm/c at LHC in quite good 
agreement also with simulations in hydrodynamical or 
kinetic transport approaches. We notice that such a val¬ 
ues correspond in one unity of rapidity to a volume of 
V - 1100 fm^ at RHIC, while at LHC V - 2500 fm^, 
which means an increase of a bit more than a factor of 
two in agreement with the estimate from pion HBT in¬ 
terferometry [^ . 

The baryon chemical potential for quarks pb = Pb/^ 
is set to reproduce the p/p ratio observed experimen¬ 
tally which means pb = 10 MeV at RHIC while at 
LHC we have approximated pb to zero. As in Ref.s 
[10 the quark and antiquark masses considered are 
'md,u = 1 X 1 ^^ = 330 MeV for light quarks, and ms = 
mg = 450 MeV for strange quarks. This implies a num¬ 
ber of quarks per unit of rapidity within the fireball at 
RHIC is Nu,d -- 230, ^ - 210 and Ns,s -- 150. At 

LHC there are 530 and Ng^g ^ 360. We also 

notice that despite the quite different total energy at both 
RHIC and LHC the energy density of the fireball (with 
radial flow energy subtracted) is about e ~ 0.7 GeV/fm^ 
in good agreement with the e nerg y density at Tc as eval¬ 
uated in lattice QCD studies [^ . 

For partons at high transverse momentum, px > 
2.5 GeV (see previous Section), we consider the minijets 
that have undergone the jet quenching mechanism. Such 
a parton distribution can be obtained from pQCD cal¬ 
culations. As in Ref. 0 we have considered the initial 
Pt distribution according to the pQCD and the thickness 
function of the Glauber model to go from pp collisions 
to AA ones.Then we have quenched the spectra with the 
modeling as in Ref. (si to reproduce the px spectrum of 
pions as observed experimentally at px ^ 8 GeV. These 
parton spectra can be parametrized at RHIC as 


dNjet J B \ 

d?px \BEpx) 


(7) 


which is the same used in 1,0 with the values given in 
the Tabl dllll The parametrization at LHC is 


dNjet 

Ai 


^4 


d?px 

i + (s)' 

^3 ' 

!+(*)“ 

4^6 


with the values of Ai given in Tabl dllll 
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K[l/GeV'^] B[GeV] n 



g 

3.18- 

liF” 

0.5 

7.11 

u, d 

9.79- 

liF” 

0.5 

6.84 

it, d 

1.89- 


0.5 

7.59 

s 

6.51 • 

10^ 

0.5 

7.36 

s 

8.02 ■ 

liF” 

0.5 

7.57 



TABLE I: Parameters for minijet parton distributions at 
midrapidity from Au+Au at y/s — 200 AGeV 



Ai 

A 2 

As 

A 4 

As 

Aq 

9 

23.46 

4.84 

8.08 

2.78 

2.79 

2.31 

u, d 

24.68 

5.11 

8.01 

0.55 

5.65 

2.56 

u, d 

23.12 

5.05 

8.21 

0.57 

5.62 

2.58 

s 

24.14 

5.11 

8.01 

0.55 

5.65 

2.56 

s 

23.12 

5.00 

8.31 

0.57 

5.62 

2.61 


TABLE II: Parameters for minijet parton distributions at 
midrapidity from Pb-Pb at y/s = 2.76 ATeV 


IV. HADRON SPECTRA AND RATIOS AT 
RHIC 

In this section, we show results for the transverse 
momentum spectra of pions, protons,antiprotons, kaons 
and Lambdas using the model described in previous sec¬ 
tions. For the coalescence contribution, we first take 
into account the effects due to gluons in the quark-gluon 
plasma by converting them to quarks and anti-quark 
pairs with probabilities according to the flavor compo- 
sitions in the quark-gluon plasma, as assumed [111 Il6l | . 
For both mesons and baryons, we include not only co¬ 
alescence of hard and soft partons as in Ref. 0 but 
also that among soft hard partons which are relevant for 
hadrons at pt ~ 3 — 4 GeV. Furthermore, we include 
both stable hadrons such as pion, proton (anti-proton), 
and kaon (anti-kaon) as well as the first excited reso¬ 
nances. With respect to 0 we have added some more 
resonance which has allowed also to verify that the inclu¬ 
sion of resonances improves the description especially at 
low Pt but does not affect significantly the intermediate 
Pt and the baryon/meson ration around the peak. 

For the pion 7r(/ = 1, J = 0) spectrum the follow¬ 
ing resonance and their decay channels are included: K* 
(1=1, J=l/2) with K* —> Ktt; p (1=1, J=l) with p —> 
tttt; w (1=0, J=l) with w —> tititt; A (1=3/2, J=3/2) 
with A —> Ntt. In Fig. [T] it is shown the predicted 
for the coalescence plus fragmentation tt^ distribution in 
Pt for RHIC in Au-|-Au collisions at y/s = 200 GeV for 
(0 — 10%) centrality; the ex per imental data are shown 
by circles 0| and triangles (3^ . We can see a general 
good agreement for pT > 1 GeV that is about the region 
where we expect the coalescence plus fragmentation ap¬ 
proach to apply. In Figlljwe show much more details of 
the calculation. By thin solid line and dashed line it is 



FIG. 1: Pion transverse momentum spectrum at RHIC in 
Au+Au collisions at y/s = 200 AGeV, 0 — 10% centrality. 
Pion production from coalescence shown by thin solid line. 
Direct pions are shown by the dash-dotted line. Pion from 
resonance decay are dashed lines for p; dash double-dot line 
for K*; double-dashed dot line for A. Pion from mini-jet 
fragmentation are the dashed line. Sum of both hadronization 
processes shown by thick solid line). Experimental data from 
PHENIX [111 113]. 


plot the contribution from pure coalescence and fragmen¬ 
tation respectively. We can see that the contribution of 
both mechanism becomes about similar for pT — 3.5 GeV 
and we will see that for the baryon instead the coales¬ 
cence will dominate in a large pt range. We can also see 
the contribution for the decay into pions coming from 
resonances, which shows that contribution from p — ^ tttt, 
dashed line, dominates up to about pT ^ S GeV which 
is about the region where anyway the fragmentation is 
starting to take over. The contribution from K* (dashed 
double-dotted line) and A (double-dashed dot line) are 
instead quite less relevant and only contribute to some 
little improvement of the description at very low pT- Of 
course for the pions it is know more or less all the hadrons 
contribute to the feed-down, see also Ref.0, but in the 
region we are interested in the resonances included are 
sufficient to have a good description of the pion spectra 
at Pt ^ 1 GeV, see FigU 

For the proton (anti-proton) (1=1/2, J=l/2) we have 
included the A (1=3/2, J=3/2) with its decay, A —>• 
Ntt, but we have also checked that adding the next 
exited state A(1440) (1=1/2, J=l/2) production and 
decay leads to negligible contribution due to the large 
suppression factor but also to the fact that about 40% 
of the decay goes through V(1440) —>■ Atttt which 
means that the proton spectrum coming from V(1440) 
is quite steep and gives some contribution only at quite 
low Pt ^ 0.5 GeV. In Fig. [2] we show the anti-proton 
transverse momentum spectrum at RHIG including co¬ 
alescence and fragmentation by thick solid line together 
with the available experimental data (circles) from Ref.s 
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FIG. 2: Antiproton transverse momentum spectrum at RHIC 
in Au+Au collisions at ^/s = 200 AGeV, 0 — 10% central¬ 
ity. Antiproton production from coalescence thin solid line. 
Direct antiproton are shown by the violet dash-dotted line; 
from A decay are the dash double-dotted line. Antiproton 
from minijet fragmentation are the dashed line. Sum of both 
hadronization processes shown by thick solid line. Experi¬ 
mental data from PHENIX [3^ . 



FIG. 3: Kaon transverse momentum spectrum at RHIC in 
Au+Au collisions at ^/s = 200 AGeV, 0 — 10% centrality. 
Kaon production from coalescence is the thin solid. Direct 
kaons are shown by the dash-dotted line. Kaons from K* 
decay are the dash-dotted line. Kaon from mini-jet fragmen¬ 
tation is the dashed line. Sum of both hadronization processes 
is shown by thick solid line. Experimental data from PHENIX 
are shown by the circles, STAR [3^ data are shown by 
the squares. 

( 3 ^ . Again the description appears to be quite good; we 
show also the relative contribution from coalescence and 
fragmentation by thin solid line and by dashed line re¬ 
spectively. We notice that for anti-protons the two mech¬ 
anism become comparable at pr — 5 GeV which means 
that the coalescence contribution is more important for 
protons with respect to pions. 


In Fig. |3]we show also the pt distribution for the 
(1=0, J=l/2) for which we have included the AT*(892) 
(1=1, J=l/2) contribution from the decay of K* —> kn. 
The impact of the next excited state Ari(1270) is again 
negligible especially at pr ^ 1 GeV. We can see that also 
for Kaons the agreement with experimental data from 
PHENIX at low pT [s^, open circles, and STAR [s^, 
squares, is fairly good in all the range of px- By dash 
double dotted line in FiglSl we see that at low pr the 
contribution from K* decay becomes important and con¬ 
tributes to have the correct slope of the spectrum as mea¬ 
sured experimentally. One can notice as in the case of pi¬ 
ons that there is some lack of yield at pr — 4 GeV where 
the fragmentation is starting to be dominant. We antici¬ 
pate that such a systematic is observed also at LHG and 
from the ratio baryon/meson we will see that it is even 
more marked for baryons and in particular for Lambda’s. 
The Pt distribution for A(1116) (1=0, J=l/2) at RHIC 



Pt [GeV] 


FIG. 4: A transverse momentum spectrum at RHIG in 
Au+Au collisions at ^/s = 200 GeV, 0 — 10% centrality. A 
production as sum of both hadronization processes is shown 
by thick solid line; from coalescence only by thin solid line. Di¬ 
rect As are shown by the dash-dotted line. A from resonance 
decay are: (dashed double dotted line), E(1385) (double 

dashed dotted line), A(1405) (dashed line with * symbols). A 
from mini-jet fragmentation is the dashed line. Experimental 
data from STAR [i^ . 

energy in Au+Au collisions at ^ = 200 GeV for cen¬ 
tral collisions is shown in Fig. |4]by thick solid line along 
with the experimental shown by circles [i^. For the A 
there are indeed several hadronic states that have a sig¬ 
nificant contribution and we have included the following 
closest resonances with their decay channels and the per¬ 
tinent branching ratios (B.R.): I]o(1193) (1=1, J=l/2); 
So —^ Ay, A(1405) (1=0, J=l/2); A(1405) —^ Ett, 
E(1385) (1=1, J=3/2) with the two main decay channels: 
E(1385) —^ Att with B.R. = 88% and E(1385) —^ Ett 
with B.R. = 11.7%. In Fig. 0]one can also find the dif¬ 
ferent contribution of these channel, see also the caption. 
Also for the Lambda the coalescence plus fragmentation 
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model appear to be able to correctly described the experi¬ 
mental data in a wide range of px- We also find similarly 
to the anti-proton that the contribution from indepen¬ 
dent fragmentation according the AKK parametrization 
becomes dominant at pr ^ 6 GeV. 

The coalescence mechanism has had the merit to nat¬ 
urally predict a baryon/meson enhancement at inter¬ 
mediate transverse momentum, especially in the region 
Pt — 2 — 4 GeV where the pj'K^, P/'^~: ^/‘2Ks reaches 
a value of the order of unity which is a strong system¬ 
atic enhancement with respect to the one observed in pp 
collisions (isj) . 



FIG. 5: Anti-Proton to positive pion ratio at RHIC from 
Au+Au collisions at ^/s = 200 GeV. The model prediction 
is the thick solid line; by dashed line the result switching off 
the radial flow of the QGP bulk matter. PHENIX data 
are shown by circles, STAR data by triangles. 

We therefore show the p/tt”, A/Ar° in the Fig.s[5]and 
El We can see that the ratio is quite well predicted from 
its rise at low px up to the peak region and then the 
falling-down behavior. However in both cases it is clear 
that in the region of pr — 5 — 7 GeV there is a lack of 
baryon yield. This is a feature that could not be ob¬ 
served when the coalescence-1- fragmentation model was 
applied a decade ago to hadronization at RHIC because 
there were no data available for proton (anti-proton) at 
Pt 4 GeV, nonetheless it appears systematically, we 
will see it also at LHC energy, so we postpone some fur¬ 
ther comment about it at the end of the next Section. 


V. HADRON SPECTRA AND RATIOS AT LHC 

We present in this Section the results from the co¬ 
alescence plus fragmentation and the comparison with 
the experimental data for Pb-fPb collisions at y/s = 
2.76 ATeV for 0 — 10% collisions. We notice that the 
results are obtained without any change or addition of 
parameters with respect to the one at RHIC in the pre¬ 
vious Section. Indeed a main aim of the present paper is 
to show that the evolution of the baryon/meson pattern 



FIG. 6: Lambda to kaon ratio at RHIG from Au-|-Au collisions 
at = 200 AGeV. The model prediction is the solid line. 
STAR data by circles [4l| . 


from RHIC to LHC can be correctly and self-consistently 
predicted. The only change with respect to RHIC is in 
the radial flow and volume of the hadronizing fireball that 
as described above is self-consistently constrained but the 
total transverse energy and the multiplicity that at LHC 
y/s = 2.76 ATeVis about a factor 2.4 and 2.7 larger with 
respect to RHIC a/s = 200 AGeV, see Section Hl-a. 



[GeV] 


FIG. 7: Pion transverse momentum spectrum from Pb-|-Pb 
collisions at ^/s = 2.76 TeV, 0 — 10% centrality. The solid 
curve includes contributions from coalescence process. The 
thick solid line is the sum of coalescence and fragmentation 
contributions. Direct pions are shown by the dash-dotted line 
with X symbols. Pion from resonance decay are the dash- 
dotted lines with symbol for p , dashed double-dotted 
lines for K*, double-dashed dotted line for A. Pions from 
mini-jet fragmentation are the dashed line. The circles are 
the experimental data from ALICE experiment [d^ . 

In Figl7]we show the total (coalescence -I- fragmenta¬ 
tion) pion spectrum by thick solid line which is a quite 
good agreement with the experimental data on tt^ from 


















ALICE Collaboration in all the pr range, except some 
lack of yield at pt ^ 0.5 GeV due to absence of all the 
resonance decays feed-down. By thin solid and dashed 
line we show the contribution from coalescence and frag¬ 
mentation respectively. We notice that the two yields 
cross at PT — 4 GeV which is about a shift of about 1 
GeV with respect to RHIC, see Fig. [T]such a shift is due 
to the larger collective flow present at LHC that shifts to 
larger px the hadrons from coalescence . The very good 
agreement of the pr distribution at LHC already shows 
that the model is able to correctly predict the evolution of 
the absolute yield and especially its px shape correctly, in 
fact no parameter of the coalescence process, essentially 
the Wigner wave function width Ap of the hadrons, has 
been modified with respect to those used in the previous 

Slprf inn fnr RTTIL! 



FIG. 8: Proton transverse momentum spectrum from Pb-|-Pb 
collisions at yfs = 2.76 ATeV, 0 — 10% centrality. The solid 
thick line is the sum of coalescence and fragmentation pro¬ 
cess. The solid thin curve shows the contribution from coa¬ 
lescence. Protons from minijet fragmentation are the dashed 
line. The circles are the experimental data from ALICE ex¬ 
periment [i^ ]. 

In FigEl it is shown the proton spectrum at LHC 
by thick solid line and compared to the experimental 
data. The agreement also in this case is very good for 
Pt > 1 GeV up to 5 GeV that is the maximum value 
with available data. At very low px as said in Section H 
we should not expect the approach to really apply. Still 
we can notice that the coalescence over predict the yield. 
The effect was partially present also at RHIG, see Fig. 
[21 The larger discrepancy at LHG could be in agree¬ 
ment with a larger annihilation of p and p that can be 
expected to be larger with respect to RHIG and could 
explain the disagreement between the statistical model 
prediction and the experimental data. 

In FigjS] we can also see that the yield of the frag¬ 
mentation process becomes comparable to the one from 
coalescence at a pt — 6 GeV which is about a 50% larger 
with respect to the pions and also a shift of about 1.5 
GeV with respect to RHIG. This is what one would ex¬ 


pect due to the larger flow at LHC and the fact that 
barvon are more affected bv it. 



FIG. 9: Kaon transverse momentum spectrum from Pb-|-Pb 
collisions at ^/s = 2.7 ATeV, 0 — 10% centrality. The dashed 
curve includes contributions from coalescence process. The 
thick solid line is the sum of coalescence and fragmentation 
contributions. Kaons from mini-jet fragmentation are the 
dashed line. The squares are the experimental data from AL¬ 
ICE experiment [43,13, [3l • 

In Fig [21 the px distribution for is shown by thick 
solid line and again one can see the good agreement with 
the experimental data [d^ in the entire range of px ■ We 
can notice that at RHIC for both pions and kaons there 
was some lack of yield in the region where the fragmenta¬ 
tion takes over, while at LHC energy for both cases the 
agreement appears quite better. This can be expected 
because the independent fragmentation function picture 
should be better constrained at energies of the order of 
TeV. 

In Fig[Tni the experimental data for the transverse 
momentum spectrum of A is shown by circles together 
with the results from the coalescence plus fragmentation 
shown by thick solid line. The different contribution from 
excited state have been calculated and have a similar rel¬ 
ative contribution as at RHIC, see Fig[3|We can see also 
for this case the good agreement for pT > 1 GeV, but 
while for p and p the data are availalble only up to 4-5 
GeV, in this case we have the availability of data up to 9 
GeV and this allows us to see that in the px region where 
the fragmentation starts to dominate, px — 6 — 7 GeV 
there is some lack of yield. This something that we could 
only marginally spot at RHIG energy mainly through 
the A/K ratio, see Fig. [5] At both RHIG and LHC 
such a lack of yield appears where coalescence becomes 
less important therefore one can say that it seems that 
the spectrum from AKK fragmentation function appears 
too flat. This may very well be because the fragmenta¬ 
tion function for baryons in general and in particular for 
A are known to be not very well constrained. On the 
other hand we notice that the fragmentation contribu¬ 
tion has been calculated for all hadrons considered with 
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FIG. 10: A transverse momentum spectrum from Pb+Pb col¬ 
lisions at ^/s = 2.76 ATeV, 0 — 10% collisions. The thin solid 
curve includes all contributions from coalescence process.A 
from mini-jet fragmentation is the dashed line. The sum of 
coalescence and fragmentation contributions is shown by the 
thick solid line. The circles are the experimental data from 
ALICE experiment [ij [45l| . 


the same = {phad/‘2.zY and this gives a global good 
description of the spectra for tt, ill at pr > 5 GeV and 
for p,p,A for px > 8 GeV. However we mention that a 
similar problem has been recently pointed-out also in pp 
collisions in Ref.sfi^: also it seems that the KKP frag¬ 
mentation can supply a better description even if still 
not satisfying. We also notice that in pp collisions the 
fragmentation function gives an excess of yield increas¬ 
ing with pt while in AA such an effect is hidden by the 
large jet quenching that is tuned to reproduce the pion 
spectra at pr — 8 GeV. On the other it still appears 
as a too flat spectrum at 3 — 8 GeV. In AA colli¬ 
sions it is likely that studies of in-medium fragmentation 
function can find a solution [13, or it could be that 
coalescence contribution should extend to large px with 
respect to the present modeling having simple spheres in 
momentum space as Wigner function and no dynamical 
role of the interaction that could lead to an extension of 
the coalescence to pair with larger relative momentum. 

We mention that recently it has been developed a pro¬ 
cess that within the coalescence plus fragmentation ap¬ 
proach could be quite important in solving this issue (49l | . 
The idea is to describe the in-medium fragmentation as 
a quark recombination of shower partons taking into ac¬ 
count also the gluon splitting into quark pairs that re¬ 
combine. Such a mechanism seems to provide a large 
contribution for baryons in the region of pr ^ 6 GeV 


In Fig llll we compare the p/tt ratio vs px shown by 
solid line with the experimental data of the ALIGE Gol- 
laboration shown by open circles. The description is 
overall quite good with some quite limited lack of pro¬ 
ton yield at p^ ~ 6 GeV. In FiglS] it is also shown by 
dashed line the p/tt ratio if the coalescence between soft 


FIG. 11: Proton to pion ratio in Pb-|-Pb collisions at ^/s = 
2.76 ATeV. The solid line is the prediction of our model. The 
empty circles are data from ALIGE experiment in collision 
Pb+Pb at 0-5% centrality [4§] . 


partons from the QGP and a mini-jet. We can see that 
the contribution is significant for px > 3GeV. The im¬ 
pact of radial flow of the soft nartons is shown bv dashed 
( 



FIG. 12: Lambda to kaon ratio in Pb+Pb collisions at ^/s = 
2.7 ATeV. The solid line is the prediction of our model. The 
circles are data from ALICE experiment in collision Pb+Pb 
at 0-5% centrality. [iSl 


In Fig. [T^ we show the results for the A/K ratio in 
comparison with the experimental data shown by circles 
[ij. We can see generally a good overall description of 
the ratio especially in the region of the peak. Comparing 
the dashed-dot line with the thick solid line we can see 
that in the peak region a quite good agreement with the 
experimental data is reached thanks to a recombination 
of thermal soft partons with a mini-jet parton. The rel¬ 
evance of such a process is present also in the EPOS ap¬ 
proach Ref. [131 , ad it is an idea that can be traced back to 
ref. M However as discussed above a.t px ^ Q GeV there 
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is a significant lack of A yield that here in a linear scale 
appears quite large. A tendency to underestimate the 
yield is visible also in even if quite smaller thanks to 
a different fragmentation scheme with respect to AKK. 
The low ratio is only slightly low because of the large 
K yield in this pr range,see FigEl and most of the dis¬ 
agreement with the data comes from the lack of yield in 
the A's distribution from fragmentation that appears too 
soft in this pr range. In Fig[T2]we also show the behavior 
of the A/K ratio if only coalescence is considered, dashed 
line, or if only fragmention is included. 


With respect to the original paper the new data ex- 
tens to higher px, especially for the A, and this as allowed 
to spot some lack of yield for baryons in the px region 
where the coalescence contribution is dying out and frag¬ 
mentation is taking over, which means at px & GeV. 
It appears that the independent fragmentation approach 
gives too hard spectra at least up to pt 8 GeV espe¬ 
cially for baryons, similarly to what has been pointed out 
for charged hadrons in pp collisions . In AA collisions 
this result seems to point to the need of an in-medium 
fragmentation process in agreement with the first results 
of a shower recombination of quarks from gluon decays 

Elisa. 


VI. SUMMARY AND CONCLUSIONS 

In this paper, we have studied the hadronization of 
the quark-gluon plasma and mini-jet partons produced 
in relativistic heavy ion collisions in terms of the parton 
coalescence plus fragmentation model. The px distri¬ 
butions of partons in the quark-gluon plasma is taken 
to have an exponential form with a temperature similar 
to the phase transition temperature, ~ 160 MeV but 
boosted by a radial flow [3. The volume and radial flow of 
the hadronizing QGP are constrained by the total mul¬ 
tiplicity and total transverse energy. Partons at higher 
Px — ISTc — 2.5 GeV are taken to be the minijets with 
power-law spectra that have undergone the jet-quenching 
process. The approach is in its core the one developed at 
RHIC energy a decade ago 0 

The aim has been to investigate if the approach already 
successful a decade ago with the first RHIG data was able 
to correctly predict the evolution of the transverse mo¬ 
mentum spectra from RHIC to LHC energy. More specif¬ 
ically the new data also at RHIC have been available in a 
wider px range especially for the baryon-to-meson ratio 
and at LHC it has become possible to verify the approach 
in a wider px range up to about 10 GeV also for the A and 
for both the p/ir and A/K ratio. We found a quite good 
agreement of both the px spectra for the main hadrons 
like TT, K,p,p, A and the baryon-to-meson ratio p/tt, A/K 
in a wide region of px up to about 10 GeV. The yield and 
slope of the spectra as well as for the height and px posi¬ 
tions of the peak in the baryon/meson ration is correctly 
predicted. This is achieved without any adjustment of 
the coalescence parameter Ap for mesons and baryons 
and hence can be considered as a real prediction of the 
model. We also note that underlying such agreement a 
key ingredient is the radial flow of the QGP matter. A 
coalescence model without radial flow could not account 
for properly neither for the baryon over meson enhance¬ 
ment nor for the single px hadron spectra. 


We remind that the coalescence has played an impor¬ 
tant role in the discussion about the QGP properties not 
only thanks to the explanation of the baryon/meson en¬ 
hancement but also because of the approximate quark 
number scaling observed at RHIG and suggested by a 
simplified approach to a pure coalescence mechanism |I4| . 
A more realistic approach to coalescence in three dimen¬ 
sion with radial flow correlations, finite hadronic wave 
function widths and resonance decays shows that about 
a 10% breaking has to be ex pec ted at intermediate pr, 
and quite larger one at low pt correctly predicting 

the experimental observation [^. At LHC energy exper¬ 
imental data show a larger breaking of the scaling with 
respect to the one observed at RHIC or predicted by more 
realistic coalescence models. However we note that the 
data discussed now are based on event-by-event analysis 
that shows the presence of higher harmonics like V 3 ,V 4 , V 5 
which also have a quite large variance. This can be ex¬ 
pected to further break the naive quark number scaling of 
the V 2 {px) but a quantitative approach requires an exten¬ 
sion of the present approach to an event-by-event Monte 
Carlo one. Such an approach is currently under develop¬ 
ment and it requires a significant longer computational 
effort to have a good statistical mapping of the various 
harmonics Vnipx) in each event and in a wide range of 
Px- However the present work, showing that also at LHC 
the main hadronic px spectra and baryon-to-meson ra¬ 
tios are well predicted, provides the motivation to extend 
the study to the anisotropic flows. 
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